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Introduction 

We  shall  consider  the  propagation  of  acoustic  disturbances  generated 
be  a  moving  point  source,  in  the  presence  of  a  plane  refracting  surface.  We 
specify  that  the  source  moves,  with  a  constant  velocity  V,  parallel  to  a 
plane  surface  separating  two  media  of  different  densities  and  sound  velocities. 
In  a  previous  paper,  the  author  (I960)  gave  a  solution  which  was  valid  when 
the  velocity  V  is  greater  than  the  sound  velocity  in  both  media,  for  the 
source  moving  along  the  interface;  this  solution  was  found  by  assuming 
conical  motion  and  then  specifying  the  nature  of  the  source.  The  method 
is  suitable  neither  in  the  case  when  the  source  is  off  the  interface,  nor 
when  its  velocity  is  subsonic  with  respect  to  one  of  the  media. 

In  the  half-space  y  >  0  we  find  a  medium  1  of  density  and  of 

sound  velocity  c  ^  ;  in  y  <  0  we  have  a  medium  2  with  density  (  P2  =  K  Pj) 


Sponsored  by  the  University  of  Melbourne,  Melbourne,  Australia  and  the 
Mathematics  Research  Center,  U.  S.  Army,  Madison,  Wisconsin,  under 
Contract  No.  DA-1 1-022-ORD-2059. 


-2- 


#  274 


and  sound  velocity  c2  .  (We  take  c2<Cj  without  losing  generality. ) 
With  i  =  1  or  2  we  may  define  velocity  potentials  <|>  ^  for  the  appropriate 
medium,  such  that  the  velocity  vector  q  and  the  pressure  change  P  are 
given  by  the  equations 


q  =  v<t>1  ,  p  =  p^j/at  , 


while  the  «|>  satisfy  the  wave  equations 


_2.  1 


(1) 


i 

The  mathematical  problem  is  to  find  the  potentials  <j>  which  while  satisfying 
the  appropriate  wave  equation  also  represent  disturbances  set  up  by  a  steadily 
moving  source,  and  which  on  the  interface  satisfy  the  conditions  that  the 
pressure  and  normal  velocity  be  continuous.  These  conditions  are  that  when 


y  =  0 


9^  8<t>2  9<l>1  94>2 

~dy  =  ~dy~  ’  9nd  “aT"  =  K  HT 

Imagine  the  point  source  to  have  a  constant  velocity  V  in  the 
z-direction.  This  steady  motion  implies  that  the  variables  t  and  z  are 
not  independent,  but  must  be  related  by  a  relation 
t  =  t  -  zV  ; 


then  the  potentials  must  satisfy  a  reduced  wave  equation 


#274 


-3- 


2  2  2 
8  <b  ,  9  4,  .  8  4> . 

■  i  Y.i  1  Ti 

0  2  +  2  ”  2  .  2 

8x  8y  y.  9t 

where 

Yt  =c  ,V(V2  -  c*)"1 

while  the  continuity  conditions  for  y  =  0  become 

9<t>2  94>1  9  4*  2 

8y  9y  ’  8  t  K8t 


(2) 


(3) 


It  is  clear  that  in  the  limiting  case  V  =  °°  ,  the  moving  point  source 
is  equivalent  to  an  infinite  transient  line  source.  It  is  also  clear  that  although 
y  i  represents  a  real  transverse  velocity  when  V  >  c  .  ,  it  is  imaginary  when 
V  <  c  ;  the  nature  of  the  reduced  wave  equation  2  changes  from  hyperbolic 
to  elliptic  and  the  form  of  the  acoustic  disturbance  will  also  change  from  one 
which  permits  discontinuities  across  characteristic  surfaces  to  one  which  is 
continuous  everywhere.  When  only  one  medium  is  present,  the  change  in  the 
nature  of  the  acoustic  disturbance  as  the  source  velocity  changes  from  super¬ 
sonic  to  subsonic  is  well  known  (See  e.  g. ,  Ward  1955).  For  the  two-medium 
problem  previous  work  has  been  restricted  to  the  examination  of  fixed  point- 
or  line-  sources  of  a  given  time  dependence.  As  well  as  mentioning  the 
classical  work  of  Sommerfeld  (e.  g.  1949)  we  may  cite  the  work  of  Cagniard 
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(1939)  and  Pekeris  (1956)  for  accounts  of  the  way  in  wluch  integral  transforms 
may  be  used.  For  the  fixed  point  source  of  arbitrary  time  dependence  there  is 
an  axis  of  symmetry,  for  the  uniform  infinite  line  source  of  arbitrary  time 
dependence  there  is  a  plane  of  symmetry,  but  when  the  point  source  is  moving 
with  a  finite  velocity,  there  is  no  such  property,  and  the  calculations  in  the 
two-medium  problem  using  integral  transforms  are  not  easy.  The  method 
described  below  is  suggested  as  a  useful  alternative. 
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Section  2;  The  point  source  moving  in  a  single  uniform  medium 

The  concept  of  a  point  source  of  unit  strength  moving  in  a  single 
medium  is  a  well  established  one  (See  e.g.  Ward  1955).  When  such  a 
source  moves  at  a  steady  supersonic  velocity  V,  the  associated  velocity 
potential,  which  satisfies  the  three-dimensional  wave  equation  with  sound 
velocity  c  ,  is 


^  ""  2  2  2  — 

2»(y  t  -  r  )  2 

2  2  "2 

within  the  conical  region  -yr>(x  +y)2  =r>0, 

where 

Y  =  cV(V2  -  c2)’1 

and 

t  =  t  -  zV  1 


(4) 


When  the  source  moves  at  a  steady  subsonic  velocity  V  the  potential  is 


*  -  Zl± 

^  ~  2  2  2  } 
4tt(  p.  t  +  r  )  2 


(5) 


where 

H  =cV{c2  -  V2)'* 

An  alternative  form  which  combines  both  these  results,  is  that  for 
any  fixed  value  V,  the  point  source  of  unit  strength  has  a  velocity  potential 
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given  in  the  region  x  >  0 ,  y  >  0  by  the  formula 

<t>  =  “2  Ri  /  - — Jdp - —  (6) 

4tt  -oo  (1  -  p  )2  [  Vt  -  px-  y(  1  -  p  )  2] 

When  V  <  c  we  have  a  convention  that  y  <  ip.  .  We  also  take  the  integration 
path  along  the  real  axis  of  the  complex  p- plane  to  pass  above  the  branch  point 
at  p  =  1  and  below  that  at  p  =  -1  .  For  points  in  other  quadrants  of  the 
transverse  (x,  y)  plane  the  signs  in  the  factor 

[Yt  ±px±y(l-p2)2]-1 

are  chosen  and  used  in  the  integral  so  as  to  make  the  limiting  two-dimensional 
solution  with  V  =  ®o f  Y  =  c  and  t  =  t  represent  a  disturbance  which  travels 
outwards  from  the  z-axis. 

The  integral  6  which  is  a  superposition  of  parametric  plane  solutions 
of  the  reduced  wave  equation 

2  2  2 
8  <(>  8  <j>  1  9  c}> 

2  +  2  ~  ~2  2 
8x  8y  y  &t 

obviously  represents  a  disturbance  which  travels  in  the  z-direction  with  a 
velocity  V.  To  show  its  equivalence  to  the  known  potentials  4  and  5  , 
we  consider  separately  the  cases  V  >  c  and  V  <  c  . 

We  change  to  polar  coordinates  in  the  transverse  plar.^,  so  that  we 


must  consider  the  integral 
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4*  “  2  ^  J  2-  ?  —  (7) 

4 n  -oo  ( 1  -  p  )2  (v  T  -  r[  P  cos  0  +  sin  0(  1  -  P  )2  ]  } 

in  the  range  0  <  0  <  ir/2  . 

Apart  from  the  branch  points  at  p  =  ±1  ,  the  only  singularities  of  this  integrand 
are  simple  poles  at  the  points  where  the  function 

{  Y  T  -  r[  P  cos  0  +  sin  0[  1  -  p  ]  2  ]  } 

vanishes.  When,  for  V>c,  y  is  a  real  velocity,  this  function  vanishes 
only  when  t  >  0  ,  and  then  at  the  points 


[  y  t  cos  0  ±  i  sin 


2  2 
0(Y  t 


r  )  2]  /r 


for  y  T  >  r  >  0 


and 


2  2  2 

p  =  [  y  t  cos  0  ±  sin  0  (  r  -  y  t  ) 


]/r 


for  r  >  y  t  >  0 


For  a  fixed  value  of  0  the  locus  of  the  complex  poles  for  all  positive  values 
of  the  ratio  x/r  is  the  branch  of  a  hyperbola  given  in  the  right-hand  half  p-plane 
by  the  equation  p  =  cosh(v  +  i0)  for  real  values  of  v.  When  r  /  0,  the 
integrand  is  0(p  2)  as  I  p| -*•<»;  the  integration  path  can  therefore  be 
shifted  by  a  rotation  in  the  positive  sense  into  the  hyperbolic  path  mentioned 


above  and  shown  in  Figure  1  without  change  of  value.  Small  deformations  must 
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be  made  in  this  path  to  avoid  the  two  conjugate  poles  in  the  sense  shown 
(for  a  given  pair  of  values  r  and  t  ),  and  the  single  real  pole  to  the 
left  of  the  hyperbola  must  also  be  noted.  Only  the  residues  at  the  two 
poles  on  the  hyperbola  contribute  to  the  integral  7  ,  and  the  formula  4 
is  easily  recovered.  The  residue  at  the  pole  on  the  real  axis  is  real  and 
makes  no  contribution  to  the  integral,  while  the  principal  value  of  the 
integral  along  the  hyperbola  is  easily  seen  to  be  imaginary  and  hence  of 
no  account.  (By  considering  separately  the  two  halves  of  the  hyperbolic 
path  we  arrive  at  an  integral  which  is  the  difference  of  complex  conjugates.  ) 

In  the  absence  of  poles  when  t  <  0,  the  integration  path  is  shifted  into  a 
horizontal  loop  between  the  points  p  =  1  and  p  =  °°  ;  there  is  no  contribution 
to  <j>  in  this  case. 

V\/hen  y  ( =  ip.)  is  imaginary  the  simple  poles  of  the  integrand  are 
for  a  given  0  at  the  points 

2  2  2  — 

p  =  [  ip.  t  cos  0±sin0(r  +p-r)2]/r  , 

these  points  being  either  in  the  upper  half  p- plane  when  t  >  0  or  in  the  lower 
half  p-plane  when  t  <  0.  The  pole  in  the  right-hand  half  plane  lies  on  the 
hyperbola  p  =  i  sinh(v  -  i0)  for  real  values  of  v,  and  it  is  to  this  hyperbola 
that  we  shift  the  integration  path.  Two  distinct  situations  arise.  When  t  >  0 
only  the  pole  in  the  first  quadrant  need  be  considered  in  the  integration.  When 
t  <  0 ,  both  the  pole  in  the  third  quadrant  and  that  in  the  fourth  quadrant  make  a 
contribution.  These  points  and  the  sense  in  which  the  poles  are  encircled  are 
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shown  in  Figures  2a  and  2b .  Again  the  principal  value  of  the  integral  along 
the  hyperbola  is  imaginary  leaving  only  the  residues  at  the  poles  to  fix  the 
potential,  and,  regardless  of  the  sign  of  t  the  result  is  given  by  equation  (5). 

To  calculate  the  strength  of  the  point  source  from  the  potential  7  we 
shall  use  an  indirect  method.  A  direct  method  involves  the  calculation  of 
volume  flux  across  a  closed  surface  which  contains  the  source.  With  the 
usual  aim  of  reducing  the  subsequent  integration  to  one  over  the  range  of 
a  single  variable  this  surface  is  chosen  differently  according  to  the  specific 
value  of  V.  When  V  <  c  an  ellipsoid  is  taken,  and  when  V  >  c  ,  a 
hyperboloid.  The  indirect  approach  avoids  the  need  for  individual  treatment. 

Starting  from  the  definition  7  for  <j>  in  the  region  r  >  0,  0  <  0  <  tt/2 
we  make  several  statements.  These  are 

1)  The  quantity  t<j,  is  a  function  only  of  the  variables  s  (  =  x/r) 


and  0  . 


2)  The  quantity  Q  which  satisfies  the  same  transverse  wave 
equation  as  <|>  and  also  the  relation  <4>  =  9Q/9t  is  given  by  the  equation 


Q  =  - 


4ir* 


R1  / 


oo  2  ^ 

I  n{y  t  -  r[  p  cos  0  +  sin  0(  1  -  p  )  2}  dp 


-oo 


7  i 

(1-pV 


It  follows  that 


t<(>  =  t9Q/9t  =  -  s9Q/9s  =  -r9Q/9r  . 
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3)  Q  is  an  integral  superposition  of  logarithmic  wave  potentials, 
and  is  therefore  the  velocity  potential  for  some  distribution  of  point-sources 
along  the  z-axis. 

4)  The  strength  when  t/0  of  such  a  line  distribution  is  found 
by  calculating  the  volume  flux  per  unit  length  S_  across  a  cylindrical 

X 

surface  which  contains  the  z-axis.  Thus 


Q 


lim  ( rdQ/di)  d6  =  -4 

r-0 


0 


lim  (  t<|>)  d0  . 
r-0 


(8) 


5)  Because  of  the  relation  between  cj>  and  Q  the  quantity 
=  8Sq/8t  is  the  volume  flux  per  unit  length  created  by  the  source 
whose  potential  is  $  .  The  total  volume  flux  created  by  this  point 
source,  found  by  integrating  with  respect  to  t  must  be  equal  to 
the  discontinuity  Sq{  t  =  0+)  -  Sg(  T  =  0_)  . 


6)  It  is  simple  to  find  the  quantity  lim  ( r<|>)  from  the  equations 

r-0 

(4)  and  (5)  .  From  (4)  and  (8)  it  follows  that  S_  =  U(t)  when 

X 

V>c  [U(t)  =1  when  t  >  0,  U(t)  =0  when  r<0],  while  from  (5) 
and  (8)  Sq  =  sgn(r)/2  [sgn(T)  =1  if  t>0,  sgn(-r)  = -1  if  t<0  J 
when  V  <  c  .  Thus  for  all  values  of  V  the  potential  7  is  associated  with 
a  source  of  unit  strength.  The  distinction  between  the  subsonic  and  super¬ 
sonic  interpretations  of  the  potential  Q  is  rather  curious.  When  V>c, 

Q  is  the  potential  of  a  semi-infinite  line  source  moving  lengthways  with 
velocity  V.  When  V  <  c  ,  Q  is  the  potential  of  an  infinite  line  singularity 
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moving  lengthways  with  velocity  V,  the  singularity  behind  the  point  z  =  Vt 
being  a  uniform  line  source,  and  that  ahead  of  this  point  being  a  uniform  line 
sink. 


The  potential  of  a  lengthways  moving  uniform  line  source  of  finite 
length  may  easily  be  found  by  simple  superposition.  Even  though  the 
fundamental  line  potentials  Q  have  different  interpretations  according 
to  the  value  of  V,  there  is  no  such  difference  in  the  case  of  finite  length. 
Thus,  for  a  length  L  of  a  line  source,  the  strength  is  defined  either  by  the 
difference  U(t)  -  U(t +LV  *)  or  the  difference  ~  sgn(T)  -  sgn  (r +LV~S 

fa  fa 

regardless  of  the  value  of  V  the  potential  Q  must  be  given  by  the  equation 

JL 


4it 


R1  /(I 


-00 


P  )  2in 


2  — 

yT  -  r[p  cos  6  +  sin  6(1  -  p  )2 

-1  2 
Y(  t  +  LV  )  -  r  [  p  cos  0  +  sin  9(  1  -  p  ) 


(9) 
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Section  3;  A  point  source  moving  parallel  to  a  plane  interface 


In  this  section  we  return  to  the  discussion  of  the  two-medium  problem. 
We  take  a  point  source  of  unit  strength  moving  in  medium  2  with  steady  velocity 
V  at  a  constant  distance  y  from  the  interface.  We  have  expressed  the 
potential  of  this  source  as  a  superposition  of  parametric  plane  waves  and 
we  are  able  to  find  the  complete  potential  by  considering  the  reflection  and 
refraction  of  these  plane  waves  and  then  writing  down  the  integral  superposition 
to  give  the  total  field. 

We  restrict  attention  to  the  region  x  >  0  .  We  have  three  distinct 


regions  of  interest  in  this  half  of  the  transverse  plane.  For  0  >  y  >  -y  , 
the  source  field  travels  in  the  positive  y-direction  and  is  reflected  in  the 
plane  y  =  0 ;  this  reflected  field  will  apparently  be  produced  at  an  image 


source  at  the  point  x  -  0 ,  y  =  -y^  and  will  travel  in  the  negative  y-direction. 
Thus  we  may  write 

(10) 


=~lt  w  / 


J  2  | 

-oo  (1  -  P  ) 


_ 1 _  _ R  (P) _ 

•  2  —  2  -■ 

-  px -  (y +  y_) (i  - p  )2  v,T-px+(y-y,)(i-p  )‘ 


when  x  >  0 ,  0  >  y  >  -y^  , 

R(  p)  being  the  reflection  coefficient  of  the  plane  wave  characterized  by  the 
parameter  p .  On  the  surface  y  =  0  both  the  primary  and  image  fields  are 

-1 

superpositions  of  plane  waves  moving  along  the  interface  with  velocity  Y2P 

2  i  _1 

and  in  which  the  variables  x  and  t  appear  in  the  function  [  y2t  -  px  -  y2(  1  -  p  )2] 
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The  parametric  plane  wave  solution  of  the  reduced  wave  equation 

2  2  2 

a  <j>1  a  <t>j  1  a  4>1 

in  medium  1  which  has  the  same  form  and  velocity  at  the  interface  and  which 
travels  in  the  positive  y-direction  must  have  the  form 

[y2t  -  px  -  y2(l-p2)1  -  y(m2  -  p2)  J]-1 
where  m  =  y2/y  .  Hence  the  potential  in  medium  1  is  given  by 
oo  Y2  T(p)  dp 

+1  - - T  W  /  - 71 - TI - ; - 7T  CU) 

4tt  -oo  (l-p  )  2(y2t  -  Px  -  y  2(  1  -  P  )2  -y(m  -pV> 

when  x>0,  y>0,  and  where  T(p)  is  the  transmission  coefficient  for  the 
parametric  plane  wave. 

For  the  potentials  10  and  11  to  satisfy  the  continuity  conditions  3 
at  the  interface,  the  elementary  plane  waves  must  also  do  so.  Hence  we 
derive  the  conditions 

«[1+  R(P)]  =T(p) 

and 

(1-P2)^[1-R(P)]  =  (m2-p2)*T(p) 


(1-pV  -  K(m2-p2)^ 

(l-p2)i  +  K(m2-p2)* 


whence 
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and 


T(P)  = 


2k  ( 1  -  p  ) 


2.  t 


(1-p2)2  +  K<m2-p2)2 


(12) 


With  these  results  only  the  field  in  the  region  x  >  0 ,  - y  >  y 
remains  to  be  given.  This  is  naturally  written  down  as  the  sum  of  the  primary 
source  field  and  of  the  image  field,  each  being  a  superposition  of  plane  waves 
which  travel  in  the  negative  y- direction. 

Thus  for  x  >  0 ,  y  >  0  , 


ra  / 


V2dP 


(13) 


2tt~  -oo  [  (1  -  p2)  2  +  K(m2  -  p2)  2  ]  [  y2  t  -  px-y2(l  -  p2)  2 -y(m2  -  p2)  2] 


for  x  >  0 ,  0  >  y  >  - y . 


«  Y2(l-P2)"  2  dp 


2.4 


2  2.4, 


and  for  x  >  0  ,  y  <  -y , 


(14) 


2 4  ,  2  2.4  ,,  2  i  .  2  2.4 

(1-p  )  +  K(m  -p  )  2  (1-p  )  2  -  K(m  -p  )  2 

I  T* 


4Tr"  -oo  [(l-p‘')2  +  K{mt‘-p‘')2]  |  Y2T-PX-(y  +  y2)(l-p2)2  Y2T-px +(y-y2>(  1-p2) 2 


1 


Y2(l-p2)’  2  dp 


2.4 


2  2.4 


(1-p  )  +K(m  -p  ) 

2  — 


(15) 

7  i  7  2- 

(1-p2)2  -  t<(m  -p  ) 2 


^2  2  J  7  -  2  2  — 

4tt  *°°  [( 1-p  )2  +  *<(m  -p  )2]  j  v2T-px+(y  +  y2)(l-p‘')2  Y2T-px+(y-y2)(  1-p 


2  - 
2--  2 


The  imputation  of  a  direction  of  motion  for  the  elementary  waves  is  strictly 
true  only  when  the  quantities  y2  and  m  are  real.  The  method  of  this  paper 
rests,  however,  on  the  basis  that  a  solution  which  is  correct  for  one  range  of 
values  of  V  is  also  correct  for  other  ranges,  since  neither  the  three  dimensional 
wave  equation  nor  the  continuity  conditions  at  the  interface  depend  on  the  value 
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of  V.  The  solutions  to  be  displayed  in  later  sections  of  this  paper  may  all, 
moreover,  be  verified  a  posteriori  . 

When  the  source  is  taken  in  medium  1  to  move  at  a  constant  height 
y  1  above  the  interface,  we  find  in  a  similar  manner  the  potentials 


V- 


oo 

Rif  - 


V2dP 


1 _ 

2  J  2—  22—  ?  ?  —  7  ( 16) 

2ir  "00  [(1-p  )2  +  k(  m  -p  )2]  [  v2T-px-y  ^ m  -p  )a  +  y(l-p  )|] 


for  x  >  0  ,  y  <  0  , 


(17) 


<b  =-  —  R1  f 
91  A  2  J 
4ir  -oo 


2  2  -  - 

oo  yz(m  -p  )  2  dp 


2  2  A 


2  — 


2  2.i 


Jk(  nT-pV+d-pV  .  K(mt'-p£')2-(l-p2)2 

i)  o  o  i  + - ^ 


2  i  A 

[(l-p  )z  +  K(m':-pV]  Y2T-px+(y-y1XmZ-p2)2  y^-px-friy^n^-p2) 


2  2 


2  2,1 


for  x  >  0  ,  0  <  y  <  y  ,  and 


♦l=-T2  “/ 


2  2  -  - 

00  Y2(m  "P  )  2  dp 


( 


2i. 


r 


4tt"  -oo  [(l-p^)2+K(m2-p2)2] 


K(m2-p2)2+(l-p2)2 


2  2  — 


(18) 

2  2  —  ?  — 

K( ni  -p  )2  -  (l-p2)2 

2  2i 


Y2T-px-(y-yi)(m  -p  )2  Y2T-px-(y+yi)(m  -p  ) 


for  x  >  0  ,  y  >  y  . 

The  six  formulae  from  equation  (13)  to  (18)  are  taken  to  define  the 
potential  field  for  any  fixed  ( real)  value  of  V.  We  shall  take  these  results 
no  further  in  this  paper;  special  attention  will,  however,  be  given  to  the 
limiting  case  with  the  point  source  moving  on  the  interface.  To  do  this  we 
can  allow  the  distances  y^  or  y2  to  approach  zero.  However  we  carry  out 
this  process,  we  have  the  expressions 
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1  r°° _ VP _ 

4*  o  ~  2  ^  t  i  ooi  2  — 

2it  -00  [(1-p  )2+K(m  -p  )2][  Y2T-Px+y(1-P  )2] 

for  x  >  0 ,  y  <  0 ,  and 


(19) 


v- 

for  x  >  0  ,  y  >  0  . 

At  this  point  it  is  worth  our  while  to  consider  the  result  of  letting  two 
sources  of  different  strengths,  one  in  each  medium,  coalesce  into  a  single 
source  moving  on  the  interface.  If  the  sum  of  the  individual  strengths  is  a 
constant  T ,  then  the  potentials  for  the  limiting  point  source  are  easily 
shown  to  be  multiples  of  the  potentials  19  and  20  by  the  factor  T.  The 
implication  is  that  the  form  of  the  point  source  potential  is  as  unique  on  the 
interface  as  it  is  in  the  interior  of  a  uniform  medium. 

We  may  add  here  the  statement  that  with  i  =  1  or  2  the  quantity 

lim  [<J>  /V]  satisfies  Laplace's  equation,  so  that  we  can  pick  out  from  the 
v-0  1 

above  equation  not  only  potentials  associated  with  the  moving  point  source, 
but  also  the  steady  potentials  associated  with  the  fixed  point  source  as  well. 
The  problem  of  the  fixed  transient  point  source  will  be  solved  elsewhere. 

In  the  remainder  of  this  paper  we  shall  examine  the  formulae  (19)  and 


00 

mf  ■ 


Y2dP 


2"“  -oo  [(l-p^) 2  +  k(  m^-p^)2]  [  Y?T-px-y(m2-p^)2  ] 


2  — 


(20) 


(  20)  in  greater  detail. 
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Section  4:  The  evaluation  of  velocity  potential  in  the  medium  with  smaller 
sound  velocity 


The  general  formula,  given  in  terms  of  polar  coordinates  r  and  0  , 

00  Y2dp 

<J>  - - ~r  R1  /  77  2  2*7  2  -i  (21) 

2 it  -oo  [(1-p  )2  +  K(m  -p  )2  (y^T-rfp  cos  6-(l-p  )  2  sin  0]  } 

represents  the  velocity  potential  in  medium  2  when  0  >  0  >  ir/Z,  but  the  process 
of  simplifying  this  integral  depends  on  which  of  fhree  possible  ranges  of  value 
the  longitudinal  velocity  has,  since  this  is  bound  to  affect  the  position  of  the 
singularities  of  the  integrand. 

When  V>  Cj>  c2  ,  the  transverse  velocities  y2  and  Yj  are  real  with 
m(=  Y2Aj)  <  1  •  The  integrand  of  (  21)  has  branch  points  at  p  =  ±1  , 

P  =  ±m  on  the  real  axis,  and  when  t  <  0  no  other  singularities.  When 
r  >  0  it  has  simple  poles  at  the  points 

2  2  2  — 

p  =  [y2tcos  0±isin0(  y2t  -r  )2]/r  when  yT>r>0 


or 

2  2  2  - 

P  =  [  y2tcos  0±sin  0(r  -  Y2T  )zjA 


when 


r  >  y2t  >  0  • 


For  a  given  value  of  0  the  locus  of  all  the  possible  complex  poles  is  the 
hyperbola  p  =  cosh  (v  -  i0)  where  v  is  real.  With  the  formal  horizontal 
integration  path  taken  to  pass  below  the  branch  points  p  =  -1  and  p  =  -m  , 
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and  above  the  points  p  =  +1  and  p  =  +m  ,  we  take  a  new  integration  path 
along  the  hyperbola,  apart  from  deformations  to  avoid  the  two  simple  conjugate 
poles  which  correspond  to  particular  values  of  r,  t  and  0  .  Then  if  we  also 
include  the  (single)  simple  pole  on  the  left  of  the  hyperbola  in  the  calculation, 
there  will  be  no  change  in  the  value  of  the  integral  when  r  £  0  since  the 
integrand  is  0(p  2)  as  lp|-»°°.  A  distinction  is  necessary,  resulting  in 
two  different  integration  paths  as  shown  in  Figures  3a  and  3b.  The  hyperbola 
in  each  case  crosses  the  real  p-axis  at  the  point  p  =  cos  0  .  When  cos  0  >  m 
the  hyperbola  passes  to  the  right  of  the  branch  point  p  =  +  m,  so  that  as  shown 
in  Figure  3b  a  horizontal  loop  integral  between  the  points  p  =  m  and  p  =  cos  0 
must  be  included  in  the  new  integration  path.  When,  as  in  Figure  3a,  cos  0  <  m, 
no  such  loop  is  needed. 

2  2  2  — 

The  pole  at  the  real  point  p  =.[y2t  cos  0  +  sin0(r  -Y2T  )2]/r  =  cos  Q  > 


say,  is  present  in  both  these  cases,  but  when  the  residue  at  this  point  is  real 
there  is  no  contribution  to  the  potential.  Thus  only  in  the  case  cos  0  >  m  and 
then  when  cos  a  >  m  can  there  be  a  contribution  to  the  potential,  because  the 
residue  now  has  an  imaginary  part. 

When  y2t  >  r  >  0  the  contribution  from  the  two  poles  on  the  hyperbola 


is 


V 


R1 


sinh  ( v-  i0) 


2  2  — 

tt(  Y“T_-r-)2  |  sinh  ( v- 19)+ K[cosh  (v-ie)-m  ]2 


2  2  2j 
2T  -r  > 

for  the  full  range  0  >  0  >  -tr/2  ,  that  is  for  both  Figure  3a  and  3b.  When 


(22) 


J  v  =  arc  cosh  y^t/ r 


the  residue  at  the  real  pole  also  is  included  we  have  the  further  contribution 
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*2=  + 


,  2  2  2  A 

"(r  -y2t  ) 


2  2  - 
sing  [ cos  a  -  m  ]2 

,2,2.  2  2 

sin  a  +  k  ( cos  a  -  m  ) 


(23) 


in  a  region  for  r  >  \  t  >  0  restricted  to  values  of  a  =  -9  +  arc.  cos  Y2T/r 
which  lie  in  the  range  1  >  cos  a  >  m  .  Equation  (22)  defines  the  potential 
in  a  sector  0  >  0  >  -  n/2  of  the  interior  of  a  cone  r  =  y2t  •  Equation  ( 23) 
defines  the  potential  in  a  'head-wave'  region  contained  between  the  cone 
r  =  y2t  >  t*le  tangent  plane  which  meets  the  interface  8=0  at  the  line 
x  =  y  jT  for  r  >  0  ,  and  the  plane  y  =  0  . 

When  for  negative  values  of  t  there  are  no  poles  of  the  integrand, 
the  horizontal  integration  path  can  be  replaced  by  a  loop  integral  round  the 
branch  points  p  =  1  and  p  =  m  ;  this  integral  is  imaginary  and  makes  no 
contribution. 

When  V<c2<c^,  the  transverse  velocities  y2  and  y^  are 
positive  imaginary, leaving  the  ratio  m  real  in  the  range  m  >  1  .  In  this 
fully  subsonic  case  the  explicit  form  of  21  to  be  evaluated  is 


1  oo  ip-2dp 

^2= - p  m  J  - 71 -  "  z  I - -  (24) 

2tr  -oo  [(1-p  )2  +  t<(m  -p  )2]  {i^2T-r[  pcos  0  -  (1-p  )2sin0]} 

with  the  formal  integration  path  along  the  real  axis  indented  so  as  to  pass 
above  the  points  p  =  1 ,  p  =  m  and  below  the  points  p  =  -1  and  p  =  -m  . 

The  integrand  24  has  simple  poles  at  the  points 
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p  =  [  i  M- 2t  cos  ®  *  s:*-n  ®(r 


2  2  i.  , 
+  Hl2T  )  l/r 


for  a  given  0  the  locus  of  poles  in  the  right-hand  half -plane  is  the  hyperbola 
p  =  i  sinh(v  +  i0)  for  real  values  of  v.  When  t  >  0  for  specified  r  and  0 
we  may  shift  the  integration  path  into  this  hyperbola  without  change  of  value, 
with  only  a  pole  in  the  first  quadrant  to  be  taken  into  account  ( See  Figure  4a). 
When  r  <  0  the  same  shift  of  path  is  possible,  but  now  poles  both  in  the 
third  and  the  fourth  quadrants  must  be  considered  (See  Figure  4b).  The  formula 
for  the  potential  turns  out  to  be  independent  of  the  sign  of  r  .  We  have  the 
result  that 


*2  = 


,  r  2  2  2  A 

2ir[r  +h2t  ]2 


R1 


cosh(v  +  i0) 


2  2  — 
cosh(v  +  i0)  +  K[m  +sinh  (v+i0)}2 


v  =  arc  sinh 


^2r  „  r°  sinh  v  cosh(v+i0)dV 

2  J  22  o  2  2  2  — 

it  0  (p.  t  -r  sinh  v)  [cosh(v  +  i0)  +K(m  +sinh  (v+i0))2]  ; 

(25) 


the  first  term  comes  from  the  residue  terms,  and  the  second  is  the  principal 
value  of  the  integral  along  the  hyperbola.  Both  terms  are  present  without 
restriction!  in  the  whole  of  the  quarter-space  r  >  0,  0  >  0  >  --rr/2  . 

In  the  intermediate  case  with  c  ^  >  V  >  c2 ,  \2  is  real,  y  ^  (  =  ip  ^) 
is  imaginary  and  the  ratio  m  (  =  -in)  is  also  imaginary.  The  branch  point 
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at  p  =  m  thus  migrates  with  change  of  V  to  the  negative  part  of  the  imaginary 
axis,  and  the  possible  integration  contour  is  also  changed.  We  are  now 
examining  the  integral 

1  r°°  v2dp 

d>  =.  —  R1  f  - 1 - 

'2,  2  *  2—22—  2  — 

2tt  *°°  [  (1  -  p  )  2-i  k(p  +n  )2]  {Y2T-r[PCOS  0-(l-p  )2sin0]}  ( 


The  position  of  the  poles  is  exactly  as  in  the  fully  supersonic  case,  but  in 
shifting  the  integration  path  to  the  hyperbola  p  =  cosh  ( v  -  10)  we  must,  in 
addition  to  considering  the  poles,  also  include  a  loop  between  the  points 
p  =  -in  and  p  =  -i«>  in  order  to  preserve  the  value  of  the  potential.  The 
full  integration  path  when  r  >  0  is  shown  in  Figure  5a  while  the  corresponding 
integration  path  when  r  <  0  when  there  are  no  poles  is  shown  in  Figure  5b. 

When  r  >  0  there  are  three  distinct  contributions.  From  the  complex 
poles  we  have  the  residue  term 


*2  =  ' 


r  2  2  2-2 

ny2  t  -r  y 


R1 


sinh(v-  i0) 


sinh(v- i0)  +  K[cosh2(v- i0)  +  n2] 2 


(27a) 


v  =  arc  cosh 


this  being  the  exact  continuation  of  equation  (  22)  within  the  cone  >  r  >  0 
The  principal  part  of  the  integral  reduces  to  an  imaginary  quantity,  so  that  we 
have  next  the  residue  at  p  =  cos  a  ,  this  being  the  head  wave  contribution 


2  2  - 
KY^  sin  a  [cos  a  +  n  ]  2 

~Tl  2  2.fr  ‘  2  '  2.  2  2- 

ir[ r  ”Y2T  J  lsln  Q  +  K  (cos  a-m  )] 


(27b) 
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with  a  =  -0  +  arccossy^/r  in  the  range  0  <  a  <  tt/2  ,  r>  y2r  <  0  .  This  head 
wave  disturbance  occupies  the  whole  region  r  >  0  outside  the  cone  r  =  V2T/* 
The  contribution  from  the  vertical  loop  is  present  for  all  positive  and 
negative  values  of  w  and  has  the  form 


k  r y2cos  0 


/ 


2  2  4 

dq(  q  -  n  ) 


r  2  2  2  2  ,  24222  2-, 

n  [1+q  *k.  (n  -q  )]  {(  y  t  +  r sin  0(  1+q  )  )  +r  q  cos  0} 


(27c) 


The  horizontal  loop  contribution  which  must  be  considered  when  t  <  0  in 
addition  that  given  b  y  equation  (  27c)  has  the  form 


k  y  r  sin  0 
6 


00 

/- 


2  4  2  2  4- 

p  -l)2  (p  +n  )2 


dp 


1  [K2(p2+h.2)+ 1-p2] [y  T-r pcos  0)  2+  (p2-l)  r2sin20] 


(  27d) 


The  general  e=xpression  for  4>2  has  now  been  simplified  in  three  distinct 
ranges  of  V.  A  featmre  which  is  held  in  common  for  all  three  ranges  is  that 
only  the  explicit  resLdue  term  is  present  on  the  axis  r  =  0  away  from  the 
source;  the  remaining  integral  terms  are  all  seen  to  vanish  on  this  axis 


away  from  the  source-. 
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Section  5:  The  evaluation  of  velocity  potential  in  the  medium  with  the 
larger  sound  velocity 


The  formula  to  be  simplified  for  various  ranges  of  the  longitudinal 
velocity  V  is,  for  0  <  0  <  it/2 


2tt 


oo 


v2dP 


4/  —  o  2  —  ( 

-oo  [(l-p^)2+t<(m  -p  )2  ]  { Y2T-r[  pcos  0  +  (m  -p  )2sin0]} 


When  V  >  c^  >  c^  the  parameters  in  the  integrand  are  real;  the  poles 
of  the  integrand  at  the  points 


P  = 


,  2  2  2  2-, 
\2t  cos  0  ±  i  sin  0  [  "y^  t  -  m  r  ] 


/r ,  for 


V  >  r  >  0  , 


lie  in  the  right-hand  half  of  the  p-plane  on  the  hyperbola  p  =  m  cosh(v  +i0) 
and  for  real  values  of  v.  This  hyperbola  crosses  the  real  axis  at  the  point 
p  =  m  cos  0  which  is  always  to  the  left  of  the  point  p  =  m  <  1  .  We  shift 
the  integration  path  to  this  hyperbola  with  indentations  at  the  poles  as  shown 
in  Figure  6.  We  note  the  presence  of  a  pole  on  the  left  of  the  hyperbola  at 
the  point 

2  2  2  2  — 

P  =  [  Y2  t  cos  6  -  sin  6  (  m  r  -  y2t  )2  ]/r 


for  r  >  Yjt  >  °  • 
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The  residue  here  is  real  and  there  is  no  contribution  to  the  potential.  The 
principal  part  of  the  integral  along  the  hyperbola  makes  no  contribution 
either,  since  the  integral  in  the  first  quadrant  combines  with  that  in  the 
fourth  quadrant  to  give  an  imaginary  quantity.  The  residue  at  the  conjugate 
poles  gives  the  result  that 


K  Y  - 


r  2  2  2,-f 

"[YjT  -r  ]2 


T  R1 


sinh  v+  i0 


2  2  _  i  ~ 

[m  cosh  (v+i0)  ]2  +  mx  sinh(v+0) 


(29) 


v  =  arc  cosh  y^/r 


for  y^t  >  r  >  0  only. 

For  r  <  0  ,  in  the  absence  of  poles  of  the  integrand,  the  integration 
path  may  be  shifted  into  a  loop  about  the  branch  points  p  =  1  and  p  =  m  ; 
this  integral  may  be  seen  to  be  imaginary,  thus  making  no  contribution  to  the 
potential. 

We  note  here  that  the  equations  (  22),  (  23)  and  (29)  determine  in 
full  the  potentials  for  a  supersonic  point  source  moving  along  the  interface; 
they  are  equal  to  these  previously  obtained  by  another  method  ( Papadopoulos 


1960). 


In  the  fully  subsonic  case  with  V  <  c  <  c  ,  and  m  =  p  /p  >  1 

1  L,  1 

we  have  the  integral 


2tr  -« 


ip2dp 


2  A 


2  2  A 


-oo  [(1-p  )2  +  K(m  -p  )2]  {i  ji^T-rf  pcos  0  +  sin  0(m  -p  )2]} 


(30) 


Here  the  poles  of  the  integrand,  at  the  points 
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Z  2  2" 

[  i|x2TCOS  G  ±  sin  0(r  2  f  h.t  )2J/r 
m 


I 

i 


lie  on  the  hyperbola  p  =  im  sinh(v  -  i 0)  for  real  values  of  v,  when 
they  are  in  the  right-hand  half  p- plane.  As  before  there  is  a  distinction 
between  the  number  of  poles  to  be  taken  into  consideration  in  the  cases 
r  >  0  and  t  <  0  ;  a  further  distinction  arises  because  when  m  sin  0  >  1  , 
the  hyperbola  p  =  im  sinh  (v  -  i© )  passes  to  the  right  of  the  branch  point 
p  =  1 .  Thus  when  m  sin  0  <  1 ,  we  shift  the  integration  path  to  the  hyperbola, 
but  when  m  sin  0  >  1,  the  same  shift  involves  the  addition  of  a  horizontal 
loop  integral  between  the  points  p  =  1  and  p  =  m  sin  0 .  The  four  cases  are 
indicated  in  Figure  7. 

It  is  found,  however,  that  the  sign  of  t  has  no  effect  on  the  formula 

for  the  potential.  Thus,  regardless  of  this  sign  and  of  the  value  of  0  in  the 

range  0  <  0  <  -rr/2  ,  we  have  the  potential 

r  n  ( 31a) 


4*1 


„  r  2  2  zA 
2u[(j.1t  +  r  ]2 


R1 


cosh(  v-  i0) 


m  k  cosh(  v-  i0)+[m  sinh(  v- i0)+ 1] 


v-  arc  sinhu^r/r 


K(x2r 


Im/ 


sinh  v 


cosh( v- i0) 


6  [  t  -r  sinh^v  ]  {itik  cosh(  v-  i© )  +[m‘'sinh‘'(  v- i0)  +  1]  2  } 


dV 


fi 

i 


i 


with  the  first  term  coming  from  the  residue  at  the  poles,  and  the  second  term 
the  principal  part  of  the  integral  along  the  hyperbola.  In  addition,  when 
m  sin  0  >  1  ,  we  have  the  contribution  from  the  horizontal  loop  integral, 
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which  reduces  to  the  term 


jx.Kr  msin6 

-VJ  - 


7  —  Z  L* 

(P  -l)z[pc°s  e  +  sine(m  -p“)2  1  dp 


2  — 


1  [p2-l  +  K2(n/-p^)]{|Jt2T2+r2[pcos  e+lm^-pV  sinG]} 


2i 


(31b) 


The  function  of  this  term  is  to  remove  a  discontinuity  in  the  quantity 
Gcj^/ae  which  appears  on  the  plane  t  =  0  in  the  first  term  of  31a 
for  values  of  0  which  make  1  >  sin  6  >  1/m  . 

Finally  in  the  intermediate  case  c  >  V  >  c  we  examine  the 

1  u 

formula 


00 


R1  /  - TT 

1  2  r/1  _2,2 


Y2dp 


--  -  »  O  i.  — 

2tt2  -oo  [(1-p14)2- iK(p2  n  )2]{y  T-r[pcos0  -isinetp'  +  nV]}  (32) 


2  — 

.  pi 


This  case  differs  from  the  earlier  ones  in  that  the  poles  of  the  integrand 
lie  at  points 


=  (y2t  cos  0 


r  2  2  2  2  , 
±isin0[r  n  +  Y7T  ]  }/r 


for  a  given  0  and  real  values  of  v  the  hyperbola  p  =  n  sinh(v  -  i0)  is 
the  locus  of  these  points,  in  the  lower  half  p-plane.  The  deformed  integration 
path  which  includes  this  hyperbola  is  shown  in  Figure  8.  Whatever  the  sign 
of  t  there  is  only  one  pole  which  makes  a  contribution  to  the  potential. 

There  are  two  other  terms  present,  the  principal  part  of  the  integral  along 
the  hyperbola,  and  the  horizontal  loop  integral.  Thus  the  potential  ^  is 


It  will  again  be  noted  that  only  the  explicit  residue  terms  are  present  on  the 
z-axis  away  from  the  point  source,  all  the  integral  terms  vanish  when  r  =  0, 
except  at  the  source  itself.  This  point  is  relevant  when  we  consider  the 
strength  of  the  source  in  this  two -medium  problem. 
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Section  6;  The  flux  associated  with  a  point  source  moving  on  a  plane 
interface 

We  return  to  the  indirect  method  for  calculating  the  source  strength 
outlined  in  Section  2.  Given  the  potential  i(>  in  the  range  0  <  0  <  tt/2  we 
define  the  quantity 


it/  2 

sn  =  -2I  lim  [  t  cj>  ]  d  0 
V1  0  r-0 


(34) 


which  represents  the  flux  per  unit  length  produced  in  medium  1  by  a  moving 
line  singularity  whose  potential  is  given  by  the  equation 

2  2  — 

oo  i  n  {y  T-r[pcos  0  +  sin  0(m  p  )2]} 

Q  “  m/  - 


2tt 


-oo 


(l-p2)2  +  K(m2-p2)2 


No  matter  what  the  value  of  V,  we  may  easily  calculate  the  value  of  S, 

V 

as  a  function  of  t  from  the  equations  (  29),  (  31),  (  33),  and  (34) 

remembering  that  only  the  explicit  residue  terms  make  a  contribution 

in  the  limit  r-*0.  In  each  case  there  is  the  same  discontinuity  in  the 

value  of  S  across  the  origin  t  =  0,  so  that  the  total  volume  flux  from 
Q1 

the  point  source  into  medium  1  is  k/(  1  +  k  ) . 

Similarly  we  have  the  quantity 


S  =  -2  j  lim  [t 4>  ]  d  0 
g2  -tt/2  r-0 


(35) 
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which  represents  the  flux  per  unit  length  produced  in  medium  2  by  the  line 
source  of  potential  Q2  given  by 

2  A 

f  n  {y2T-r[pcos  0  -  sin  6(l-p  )2]}dp 
(1-p2)2  +  K(m2-p2)2 

We  find  from  the  expressions  (22),  (25),  (27)  and  (35)  that  the  volume 
flux  created  in  medium  2  by  the  moving  point  source  is  1/(1  +  k). 

The  total  volume  flux  produced  by  the  point  source  is  unity,  it  is 
apportioned  in  a  definite  ratio  between  the  two  media  such  that  the  rate  at 
which  mass  is  created  by  the  source  is  the  same  for  each  medium  regardless 
of  the  (constant)  velocity  of  the  source. 
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Figure  1 

The  original  and  subsequent  integration  paths  for 
equation  (7)  in  the  supersonic  case  (t  >  0)  . 


Figure  2a 


Figure  2b 


The  original  and  subsequent  integration 
paths  for  equation  (7)  in  the  subsonic 
case  (  t  >  0)  . 


The  original  and  subsequent  integration 
paths  for  equation  (7)  in  the  subsonic 
case  (  t  <  0)  . 


Figure  3a 


The  original  and  subsequent  integration 
paths  for  equation  (21)  in  the  fully 
supersonic  case  (  t  >  0  ,  cos  6  <  m)  . 


Figure  4a 

The  original  and  subsequent  integration 
paths  for  equation  (  24)  in  the  fully 
subsonic  case  (t  >  0)  . 


Figure  3b 

The  subsequent  integration  path  for 
equation  (  21)  in  the  fully  supersonic 
case  (  t  >  0 ,  cos  0  >  m)  . 


Figure  4b 

The  original  and  subsequent  integration 
paths  for  equation  (  24)  in  the  fully 
subsonic  case  (t  <  0)  . 


Figure  5a 


Figure  5b 


The  original  and  subsequent  integration 
paths  for  equation  (26)  in  the  case  of 
intermediate  velocity  (t  >  0)  . 


The  subsequent  integration  path  for 
equation  (  26)  in  the  case  of  inter¬ 
mediate  velocity  (t  <  0)  . 


Figure  6 

The  initial  and  subsequent  integration  paths  for 
equation  (  28)  in  the  fully  supersonic  case  (  t  >  0)  . 
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Figure  7  a 

The  initial  and  subsequent  integration 
paths  for  equation  (30)  in  the  fully 
subsonic  case  (  t  >  0 ,  m  sin  0  <  1)  . 
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Figure  7c 

The  subsequent  integration  path  for 
equation  (30)  in  the  fully  subsonic 
case  (t  >  0,  m  sin  0  >  1)  . 


Figure  7b 

The  initial  and  subsequent  integration 
paths  for  equation  (30)  in  the  fully 
subsonic  case  ( r  <  0 ,  m  sin  0  <  1)  . 
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Figure  7d 

The  subsequent  integration  path  for 
equation  (30)  in  the  fully  subsonic 
case  (  t  <  0 ,  m  sin  0  >  1)  . 
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Figure  8 

The  subsequent  integration  path  for  equation  (  32) 
in  the  case  of  intermediate  velocity . 
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